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Abstract. It is possible to generalize the fruitful interaction between (real or complex) 
Jacobi matrices, orthogonal polynomials and Fade approximants at infinity by considering 
rational interpolants, (bi-)orthogonal rational functions and linear pencils zB — Aof two 
tridiagonal matrices A, B, following Spiridonov and Zhedanov. 

In the present paper, beside revisiting the underlying generalized Favard theorem, we 
suggest a new criterion for the resolvent set of this linear pencil in terms of the underlying 
associated rational functions. This enables us to generalize several convergence results for 
Fade approximants in terms of complex Jacobi matrices to the more general case of conver- 
gence of rational interpolants in terms of the linear pencil. We also study generalizations 
of the Darboux transformations and the link to biorthogonal rational functions. Finally, 
for a Markov function and for pairwise conjugate interpolation points tending to oo, we 
compute explicitly the spectrum and the numerical range of the underlying linear pencil. 



1. Introduction 



The connection with Jacobi matrices has led to numerous appUcations of spectral tech- 
niques for self-adjoint operators in the theory of orthogonal polynomials on the real line 
and Pade approximation. In order to give an idea of these interactions consider a Markov 
function of the form 



(p{z) 



dfi{t) 



t ' 



where a, b are real numbers and dfi{t) is a probability measure, that is, dij.{t) = 1. It 
is well known IT], ll33l that one can expand such a Markov function tp into the following 
continued fraction 
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where bj , aj E K, aj > 0. Continued fractions of the form dl.lb are called J-fractions Il23l 
l33l . To the continued fraction dl.lb one can associate a Jacobi matrix A acting in the space 
of square summable sequences and its truncation A^Q-n] 
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Then it is known that ip{z) = {{zl — A) ^cq, eg), and the nth convergent of the above 
continued fraction is given by 



((z/- ^[0:„_i]) eo,eo 



n-2 



qn{z) \z-aQ \z-an-i 

where the column vector eo = (1, 0, ... )^ is the first canonical vector of suitable size, g„ 
are orthogonal polynomials with respect to d/i, and p„ are polynomials of the second kind, 
see 111126^ 27). It is elementary fact of the continued fraction theory that 

9(^}~'r^=o(^) , (1.2, 



see for instance ||T1|4]|23|. Relation ( 11.2b means that the rational function p„/g„ is the nth 
diagonal Pade approximant to ip at infinity. Consequently, the locally uniform convergence 
of diagonal Pade approximants appears as the strong resolvent convergence of the finite 
matrix approximations A[o:„]. For instance, one knows that Pn/qn — > <y5 in capacity in 
the resolvent set of A given by the complement of the support of fi, and locally uniformly 
outside the numerical range of A given by the convex hull of the spectrum of A, see for 
instance [30|. Besides, it should be mentioned here that an operator approach for prov- 
ing convergence of Pade approximants for rational perturbations of Markov functions was 
proposed in ifTTl . see also ifTSl . 

If ip is no longer a Markov function but has distinct nth diagonal Pade approximants 
at infinity, we may still recover these approximants as convergents of a continued fraction 
of type (II. lb . but now in general aj, bj £ C, aj ^ 0, see ll33l . that is, A becomes com- 
plex symmetric, called a complex Jacobi matrix. There is no longer a natural candidate 
for the spectrum of A, but it is still possible to characterize the spectrum in terms of some 
asymptotic behavior of the Pade denominators qn{z) and the linearized error functions 
rn{z) = qn{z)(j){z) — pn{z) [S] [12] (TO) , See also ||8][T7][T6| for more general banded ma- 
trices. Convergence outside the numerical range was established in lfT2l . and convergence 
in capacity in the outer connected component of the resolvent set in ifTTl . We refer the 
reader to |[9l for some recent summary on complex Jacobi matrices, including some open 
questions partially solved in I?]. 

The goal of this paper is to generalize several of the above results to the case of multi- 
point Pade approximants. 

Definition 1.1 (f4l)- The [ni\n2] multipoint Pade approximant (or rational interpolant) 
for a function tp at the points {zk\^^i is defined as the ratio p/q of two polynomials 
p and q of degree at most ni and n2, respectively, such that (pq — p vanishes at 
Zl, Z2, z„j+„2+i counting multiplicities. 

It is easy to see that the degree and interpolation conditions lead to a homogeneous 
system of linear equations, and thus an [ni |n2] multipoint Pade approximant exists. Also, 
one may show uniqueness of the fraction p/g. However, since the denominator may vanish 
at some of the interpolation points, it may happen that the fraction p/(7 does not interpolate 
ip at some point Zk, usually referred to as an unattainable point. 

Under some regularity conditions, the [n — l|n] multipoint Pade approximants of p may 
be written as nth convergents of a continued fraction of the form 

1 I al{z ~ zi){z - Z2)\ a\{z - z^){z - zM 



Z — Oo Z — Oi Z — 02 

the odd part of a Thiele continued fraction Continued fractions of this type are re- 
ferred to as MP-fractions in 11211 and as i?//-fractions in Il22l . In particular, the authors 
study in 11211 Theorem 4.4] and l!22l Theorem 3.5] some analog of Favard's theorem and 
the link with orthogonal rational functions. Spiridonov and one of the authors 1291 l34l 
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showed that such continued fractions are related not to a single Jacobi matrix but to a 
pencil zB — A with tridiagonal matrices A, B. Various links to bi-orthogonal rational 
functions have been presented in |[34l and ifTSl . In particular, in ifTSl Theorem 6.2] the au- 
thors present an operator-theoretic proof for the Markov convergence theorem multipoint 
Pade approximants Il20ll based on spectral properties of the pencil zB — A. 

The aim of this paper is to present further convergence results for the continued fraction 
(11.3b . both in the resolvent set and outside the numerical range of the tridiagonal linear 
pencil zB — A. To be more precise, denote by l"^ ~ ^fo oo) Hilbert space of complex 
square summable sequences (xq, zi, . . . )^ with the usual inner product 

oo 

j=o 

We will restrict our attention to the case of tridiagonal matrices A, B with bounded en- 
tries, in which case we may identify via usual matrix product the matrices A and B with 
bounded operators acting in Notice that many algebraic relations remain true in the 
unbounded case as well. However, already the simpler case of bounded pencils allows to 
describe the main ideas of how to generalize results from the classical theory of orthogonal 
polynomials to the theory of biorthogonal rational functions as well as to the multipoint 
Pade approximation. 

The remainder of the paper is organized as follows; we start from a general bounded 
MP-fraction and introduce the associated linear pencils together with the rational solu- 
tions of some underlying three term recurrence relations in In m.2\ by generalizing 
previous work of Aptekarev, Kaliaguine & Van Assche 13j we show how the asymptotic 
behavior of these rational solutions allows to decide whether the linear pencil zB — ^ is 
boundedly invertible. In particular, we deduce in Corollarv l2.6l the pointwise convergence 
of at least a subsequence of our multipoint Pade approximants towards what is called the 
TO-function (or Weyl function) of the linear pencil. Subsequently, we present in Theo- 
rem l2.10l of ii2.3l an alternate proof for a Favard-type theorem based on orthogonality prop- 
erties of associated rational functions, which yields in Corollary 12.121 a simple proof for 
the fact that the convergents of our continued fractions are indeed multipoint Pade approx- 
imants of the TO-function of our linear pencil. In fj3]we generalize the above-mentioned 
results of lfT2l Theorem 3.10], ifTTl Theorem 3.1], and ifTTl Theorem 4.4], on the conver- 
gence of Pade approximants at infinity in terms of complex Jacobi matrices to the more 
general case of multi-point Pade approximants in terms of linear pencils zB — A. The 
aim of 21 is to explore LU and U L decompositions of our linear pencil, and the link to 
biorthogonal rational functions. This naturally leads us to consider generalizations of the 
Darboux transformations of fl3^. Finally, we generalize in |2]the findings described in the 
begining of this section, namely, if we start with a Markov function and pairwise conjugate 
interpolation points tending to infinity, then the spectrum of our linear pencil is still the 
support of the underlying measure, and the numerical range equals its convex hull. 

2. Continued fractions, linear pencils, and their resolvents 

In this section we show the links between continued fractions in question and linear 
pencils. Moreover, we prove a Favard type result for the corresponding recurrence relation. 

2.1. Linear pencils. Let us consider a continued fraction of the form 

1 I a^iz)a§{z)\ af(z)af(z)| 
I /3i(z) I p^iz) 

where /3n,ctn, a^' are polynomials of degree at most 1 and not identically zero. Next, 
denote by Cn{z) the nth convergent of this continued fraction obtained by taking only 
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the first n terms in (12.1b . then the well-known theory of continued fractions tells us that 
Cn{z) = Pn{z) / qn{z), where the polynomials p„ of degree <n — \ and qn of degree < n 
are obtained as solutions of the three-term recurrence relation 

Vn+i ^ I3n{z)yn - a^i-i{z)an-i{z)yn^i, n = 0,l,2,... (2.2) 
by means of the initial conditions (setting ~ = 1 for convenience) 

qo{z)^l, g-i(z) = 0, p_i(z) = -l. (2.3) 

Using ( 12.21 ) and ( 12.3b one easily verifies by recurrence that 

qn{z) = det(zB[0:„-l] - ^[0:n-l]), Pn{z) = det(2;B[i.„_i] - ^[i:„_i]). (2.4) 

By Cramer's rule, this implies the following formula for the convergents 

P7l{z) 



Cn{z) 



qn{z) 



((2:-B[0:„-l] - ^[0:n-l]) ^Gq, 60). 



(2.5) 



By induction, one also easily shows the Liouville-Ostrogradsky formula (for the classical 
case, see H] p. 9 formula (1.15)]) 



Tl-l 



Pn+i{z)qn{z) - pn{z)qn+i{z) = ]J a^{z)akiz), n = 0, 1,2, 



fc=0 

For a complex number (f>{z), the sequence defined by 

r„(z) := 0(z)(7„(z) - p„{z), 
gives another solution of ( 12.2b with initial conditions 

ro{z)^4){z), r_i(z)=:l. 



(2.6) 



(2.7) 



(2.8) 



We will refer to r„ as linearized error (or function of the second kind) since, from the 
Pincherle Theorem 1231 Theorem 5.7], the continued fraction ( 12.1b has a limit 0(z) iff 
r„(z) is a minimial solution of the recurrence relation (12.2b . 

It will be convenient to write the polynomials a^', a^, and (3j occurring in (12.1b in the 
form of the tridiagonal infinite Unear pencil 

f Mz) -a^iz) .. 

-a^{z) /3i(z) -«f(z) 

-a^iz) hiz) a§{z) 



zB-A = 



\ 



(2.9) 



7 



with the two tridiagonal infinite matrices A = {ai,j)f°,Q and B = (6 



For a J- 



fraction, we obtain the linear pencil z — A with a tridiagonal matrix A ||T] (see also ||9l). 
In the case of J-fractions it is also known that we may write the eigenvalue equation 
Ay = zy for some infinite column vector y in terms of normalized counterparts of the 
monic polynomials qn{z) (namely the corresponding orthonormal OP). Notice that the 
product Ay is defined for y not necessarily an element of since for each component 
there are only a finite number of non-zero terms. For the linear pencil zB — A we can 
analogously write the similar eigenvalue equations 



Aq^iz) ^ zBq^iz), q^{z)A^ zq^{z)B, 



(2.10) 



with an infinite column vector q^{z) = ((/^(z), qf'{z), . . . )^ and an infinite row vector 
q^(z) ~ (qq (z), (z), . . . ). Here g^(z) and q^'iz) are rational functions obtained from 
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g„(z) by scaling with a product of linear polynomials. Indeed, defining Qniz), Pn{z), and 
r^{z) via 



quiz) = 



Ji-1 



Pn{z) 



n "^(^) n "^(^ 

fc=0 fc=0 

leads us to three solutions of the recurrence relation 
In the similar way, we see that 

qn{z) 



n "^(^) 



(2.11) 



fc=0 



0, 



0,1,2, 



Pn{z) 



fc=0 



rn{z) 
_-l ■ 



(2.12) 



(2.13) 



fc=0 



0, 



0,1,2,. 



(2.14) 



are three solutions of the recurrence relation 

a^(^)yn+i - Pn{z)y^ + a^_^{z)y^_^ 

Now, it is immediate to see by taking into account the initial conditions ( I2.3l l that the 
identities ( I2.12l i and ( 12.14b reduce to the formal spectral equations (12.10b . 
It should be also noted that we formally have 

A) = -£,([, {zB - A)p"{z) = ~eo. (2.15) 



Remark 2.1. There are many degrees of freedom in going from a continued fraction 
(12.1b to a linear pencil zB — A. For instance, for the special case of dcg /3„ = 1 and 
dcg ^ ~ deg a^- for all n > 0, the above approach leads a priori to diagonal B and 
tridiagonal A without any further symmetry properties. However, by applying an equiva- 
lence transformation to ( 12.1b we can make the polynomials /3„ monic, implying that B is 
the identity matrix. Moreover, we can choose = a^, i.e., A becomes complex sym- 
metric (also called a complex Jacobi matrix). In this case, qf^ = are known to be the 
corresponding formal orthonormal polynomials, whereas g„ is the associated monic coun- 
terpart. We will return to this scaling and normalization freedom in the last section. □ 

2.2. m-functions of linear pencils and the resolvent. In accordance with the Jacobi case 
of B being the identity, we define the resolvent set p{A, B) of the linear pencil zB — Ato 
be the set of z G C such that zB — A has a bounded inverse. The following simple example 
shows that p{A, B) can be of arbitrary shape. 



Example 2.2. For 29,2:1, 
Diz — D2 = dia 



G C, consider the diagonal pencil 



1 + \^n\J „=0,1,2, 

together with the tridiagonal pencil 



/ 1 -1/2 
1 -1/2 



Bz- A = U*{Diz - D2)U, 



U 











1 





-1/2 



V 



with bounded A, B. Then according to ( IZ9T l. (irTTT i, and (imi l. 

a^(z) = a^{z) = \{^y Quiz) = = 2". 
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Since U is boundedly invertible, we find that p{A, B) = p{D2, Di) = C \ S with S the 
closure of {zo, zi, . . .}, which could be any closed set of the complex plane. In particular, 
p{A, B) can be empty. □ 



Notice that, if in addition B is boundly invertible (which for instance in Example 12.21 
is not necessarily true) then p{A, B) — p{B^^A) — p{AB~^). However, in general we 
loose for the last two operators the link with tridiagonal matrices and three-term recurren- 
cies, and thus we prefer to argue in terms of pencils. 

Aptekarev et al. |3j Theorem 1] showed that a bounded tridiagonal matrix has a bounded 
inverse if and only if the above solutions of the recurrencies ( 12.12b and ( 12.14b have a par- 
ticular asymptotic behavior. In our setting, their findings (see also the slight improvement 
given in fT2l Theorem 2.1]) read as follows. 



Theorem 2.3 (ID). Suppose that A, B are bounded, and consider for z G C the matrix 
R{z) with entries 

r rf{z)q^[z) = {qfiz)cbiz)-pf{zMiz) if J > k, 
- 1 (z)r^(z) = gf(z)(g,^(z)0(z) ifj < k. 

Then z £ p{A, B) if and only if there exists (t>{z) G C and constants 7(2;) > 0, 5{z) G 
(0,1) such that 

<7W'5W''"", j,fc = 0,l,.... (2.16) 
In this case, R{z)j,k = {{zB — A)^^ek, ej), in particular, 4>{z) is uniquely given by 
= i?(z)o,o = {{zB - Ay^eo, eo). 

For bounded complex Jacobi matrices (B = I and q^{z) ~ q^{z)), it was shown in 
ifTOl that z G p{A, B) can be characterized only in terms of the asymptotic behavior of 
the denominators q^i^z) = q^iz)- As we see from Example 12. 21 this is no longer true for 
bounded tridiagonal pencils. 

For the sake of completeness, we will give below the ideas of the proof of Theorem l2.3l 
Let us first discuss some immediate consequences. 

Remark 2.4. For the particular case of Jacobi matrices (that \s B — I), the above formulas 
for the entries of the resolvent, also referred to as Green's functions, have been known for 
a long time, see for instance the recent book Il28l Section 4.4]. Our linear pencil formalism 
also includes so-called CMV matrices occurring in the study of orthogonal polynomials 
on the unit circle, see 1281 Section 4.2], here A, B are not only tridiagonal but in addition 
block-diagonal, with unitary blocks. Again, the formulas for the Green's functions given 
in 1281 are a special case of Theorem l2.3l We also refer the reader to Q for recent findings 
for the special case of multipoint Schur functions : here the roots of and are related 
through reflexion across the unit circle. □ 

A basic object in Theorem l2.3l and in the rest of the paper is following. 

Definition 2.5. The function 

m(z) = ((zB-^)-ieo,eo), zep{A,B) (2.17) 
will be called the m-f unction (or Weyl function) of the linear pencil zB — A. 



Comparing with ( 12.51 ) we are left with the central question whether the m-function 

Pn{z)/qn{z) = ((2S[0:„_i] - ^[0:n- 1] ) " ^ Cq , Gq) of the finite pCUCil 2S[0:„_i] - A[o,n-l] 

converges for n — > oo to the m-function of the infinte pencil zB — A. 

We learn from Theorem l2.3l that the linearized errors (z) = R{z)o^n = {z)'m{z) — 
Pn{z) and r^{z) ~ R(z)n.o = Qn{z)m{z) — p^{z) tend to zero with a geometric rate 

limsup|r^(z)|i/" < 1, limsup|r^(z)|i/" < 1, z G (2.18) 
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Following exactly the lines of Aptekarev et al. 13] Theorem 2 and Corollary 3] we obtain 
the following result on point-wise convergence of a subequence. 



Pniz) 



1/n 



< 1. 



Corollary 2.6. We have for z e p{A, B) 

limsup |g,^(2:)|^/" > 1, limsup > 1, liminf m{z) 

n— ^oo n— )-oo n— »oo Q^yZ) 

Proof. Using (12. lib and (12.131 1. the Liouville-Ostrogradsky formula (12.6b takes the follow- 
ing form 

«n(^)9^+l(^)^f (^) - «n(^)'Z^(^)^^+l(^) = 1- (2.19) 

Since sup„ max{|a^(z)|, |a^(z)|} < oo by assumption on A, B, relation (12.18b together 
with the Cauchy-Schwarz inequality implies that 

liminf[|g,^(z)P + |g,^^^(z)p]V(2n)>i, 

n— ^oo 

implying our first claim. The second is established using similar techniques, and the third 
by writing m(z) -p„(z)/g„(z) = r!^{z)/q^^{z). □ 

By having a closer look at the proof, we see that we have pointwise convergence for a 
quite dense subsequence, namely for p„-|_e„ /9n+e„ for n > with suitable e„ G {0, 1}. 
We will show in Theorem l3 . 5 I below that this point- wise convergence result can be replaced 
by a uniform convergence result in neighborhoods of an element of p{A, B). 

In the remainder of this subsection we present the main lines of the proof of Theo- 
rem l2.3l The first step consists in showing that our infinite matrix R{z) is a formal left and 
right inverse for zB — A, compare with |[33l Section 60 and Section 61] for the case of 
complex Jacobi matrices. 

Lemma 2.7. For any value of (f>{z), the formal matrix products R{z)(zB — A) and (zB — 
A)R{z) give the identity matrix. 

Proof. We will concentrate on the first identity, the second is following along the same 
lines . Write shorter 

110:3] - iqo,---,qf,o,o,---) 

and similarly p^.^j and r^.^j for the row vectors built with the other solutions of the recur- 
rence ( 12.14b . Then 

pfo.,,^iz)izB -A) = -e^ + afiz)pf^,{z)ej + afpfiz)ej^,, (2.20) 

q^O:j]iz){zB ^A)^ af{z)q!f^,{z)ej + af q^{z)ej+,. (2.21) 
In view of ( 12.6b . ( 12.1 lb . and ( 12.131 ), one obtains 

(gf (z)pf„^^.](z) -pf (^)gfo^^.j(z))(zB - A) - ej - q;f{z)ej. 
In addition, from ( 12.10b and ( 12.15b we have that 

{q'^{z)<j>{z)-p^{z)){zB-A)^eJ. 
A combination of the last two equations shows that, for all j > 0, 

(i?(2),- 0, i?(2),,2, . . . ){zB -A) = eJ, 

as claimed above. □ 
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Proof of Theorem \2.3\ Let z e p{A, B). Then, according to Lemma IZTl R{z) is indeed 
the matrix representation of the bounded operator [zB—A)^^. We get the decay rate ( 12.16b 
of the entries of R{z) from ifTSl Theorem 2.4] using the fact that R{z) is the inverse of a 
bounded tridiagonal matrix. 

Suppose now that 0(z) e C is such that ( 12.16b is satisfied. Then, using the same argu- 
ments as in Q we have that R{z) represents a bounded operator in which by Lemma lZTl 
is a left and right inverse of zB — A. Hence z e p{A^ B). □ 

Remark 2.8. The essential tool in the proof of Theorem 12.31 was the decay rate ( 12.16b 
of entries of the inverse of a bounded tridiagonal matrix. In order to specify the rate of 
convergence, for instance in Corollarv l2.6l it is interesting to quote from ifTSl Theorem 2.4] 
possible values of 7(2), 5{z) in terms of the condition number 

k[z) = \\zB-A\\ \\{zB-A)-^\\ > 1 

being obviously continuous in p{A, B), compare with lfT2l Lemma 3.3], 

j n{z)-l 2,\\{zB-A)-^\\ r {I + K{z)f ^ 
^(^) = V^<^' 2.[z) I- 

□ 

2.3. Biorthogonal rational functions and a Favard theorem. Our explicit formulas for 
the entries of the resolvent allow for a simple proof of biorthogonality for the denominators 
and q^, and in addition an explicit formula for the linear functional of orthogonality 
discussed by Ismail and Masson L22J . This generalizes the classical case of B = I and a 
selfadjoint Jacobi matrix A (l\ where it is well-known that, for j ^ k, 

{qj{A)eo,qk{A)eo) = 0. 

As a consequence, we obtain a simple proof of the fact that the 71th convergent of ( 12.1b is 
indeed an [n — 1 |n]th multipoint Pade approximant of the m-function. 

In this subsection we denote for fc = 0, 1, 2, ... by Z2k+i (and by Z2fe+2) the root of 
(and of a^, respectively), where we put Z2fc+i = 00 (and Z2fc+2 — 00) if (and a^) is 
of degree 0. Similar to ll22l we suppose for convenience that zi, Z2, ... S p{A, B). More 
precisely, we suppose that there exists a domain T^xt with compact boundary forming a 
Jordan curve such that 

zi, Z2, ... e T,,t C Clos(re,t) C p{A, B), (2.22) 

where Clos( ) denotes the closure. The case = 00 needs special care: notice that 
00 G p{A, B) if and only if B has a bounded inverse, in which case we will also suppose 
that 00 e Text- The boundary F of Text is orientated such that Text is on the right of F, 
implying that 

( ^ W 5(c) 

for z £ Text and any function g being analytic in p{A, B) and, if 00 G p{A, B), vanishing 
at infinity. 

We start by establishing an integral formula for the entries of the resolvent. 
Lemma 2.9. Under the assumption (12.221 ). we have for z G Text and j, fc = 0, 1, 2, ... 

R{z),.u ^ {{zB-A)-^eu,e,) ^ ±- J^qf {C)q^{C)^dC 
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Proof. We will consider only the case j > k, the case j < k is similar. Both the resolvent 
and Rj^k are analytic in p{A, B), and vanishing at infinity provided that oo E p{A, B). 
Using the explicit formula for R{z)j^k derived in Theorem l2.3l we get for z € T^xt 



It remains to show that the last integral equals zero. Denote by a connected component 
of C \ Clos(rea;t). If is bouudcd, then, by assumption ( 12.221 ). all poles of the rational 
function ^ i— pfiO^kiC) /{^ ~ ^6 outside of Clos(O), and hence the integral over dfl 
is zero. If 51 is unbounded, then by the above assumption on Text we may conclude that 
Qo ^ p{A, B), implying that all ze are finite. It follows from ( 12.1 lb and ( 12.131 1 that all poles 
of the rational function C vfiOlk (C)/(^ ~ C) ^6 outside of Clos(fi), and this function 
does vanish at oo. Hence again the integral over is zero. □ 

We are now prepared to state and to give a new constructive proof of the Favard type 
Theorems ll22l Theorem 2. 1 and Theorem 3.5] of Ismail and Masson. 

Theorem 2.10. Under the assumption (12.221 1. define for g e C{T) the linear functional 

6(.9) = ^^5(CMC)dC, 

then we have the following biorthogonality relations: for any n > 1 and for any polyno- 
mial p of degree < n there holds 

Proof. We again only show the first relation, the second follows by symmetry. Observe 
first that af^_i{z2n-i) = implies that zB — A is upper block-diagonal. Since Z2ri-i G 
p{A, B) by (I2.22I 1. we obtain for the resolvent {z2n-iB — A)^^ the block matrix repre- 
sentation 



(Z2,i-li?[0:„-l] — ^[0:n-ll) ^ 


* 








In particular, comparing with (12.4b it follows that qk{z2k-i) 7^ for fc = 0, 1, n — 1 (or 
degqk = k provided that Z2k-i = 00), and 

R{z2n~l)n,k = 0, fc = 0, 1, ...,71 - 1 

(or lim^^oo zR{z)n,k = in the case Z2„-i = 00). The first relation impUes that 
span < : fc = 0, 1, n - 1 > = l — z — ■ 'i^SP < ^ 



and the second combined with Lemma lZ9] that 

®te^)=0' fc = 0,l,...,n-l, 
as claimed in Theorem l2.10l □ 

Remark 2.11. In the statement of Theorem 12. 101 one recovers the jn-function as a gener- 
ating function for the linear functional of orthogonality, since 

z ^ <5q(^—^-^ ^ m{z), zeVext- 
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Suppose in addition that oo G p{A,B), and thus B has a bounded inverse. Then Cauchy's 
theorem gives the normahsation ©(1) = m'(oo) = {B^^cq, cq), and for £ > 

Similarly, for Zk E Text and ^ > we have that 



Using a partial fraction decomposition, we obtain for any polynomial p (of degree < 2?! if 
oo ^ p{A, B)) the even simpler formula 

e{r)^{B-'riAB-')eo,eo), ~ ^ 



□ 

The orthogonality relations of Theorem 12. 1 01 allow now to show in a simple way that 
the convergents of our continued fraction ( 12.11 ) are indeed multipoint Pade approximants. 

Corollary 2.12. Under the assumption ( I2.22l i, for any n > 0, the rational function Pn/ q_n 
is an [n — 1 |n] multipoint Fade approximant of the m-function of the pencil zB — Aat the 
points zi, Z2n counting multiplicities. 

Proof. Relation ( 12.4b shows that p„, are polynomials of degree at most n — 1, and n, 
respectively, and from the proof of Theorem l2.10l we know that g„(z2n-i) 7^ 0, hence g„ 
is non-trivial. 

The interpolation conditions for a Cauchy transform (or more generally for a generating 
function of a linear functional) are known to translate to orthogonality relations with vary- 
ing weights, see for instance ll30l Lemma 6.1.2]. Since = {mqn — Pn)/{oiQ ■ ■ . a^-i) 
is analytic in Text (and vanishes at 00 if 00 € Text), we only have to show that ui := 
r^/{aQ . . . ct^^i) is analytic in Text, and, provided that 00 G Text, its expansion at 00 
starts with a term z~"~^. 

Denote by zi, the finite points out of zi, 2:3, Z2k-i with fc < n. If £{k) > 1, 
define 

a^(z)...a,^_i(z) 
P[z) = 7 ^ 7 r, 

[Z - Z£(l)) . . . (Z - Z£(fe)) 

a polynomial of degree < n. Arguing as in Lemma 12.91 and using the Hermite integral 
formulas for divided differences we find that 

where in the last step we have applied the orthogonality relation of Theorem l2.10l Hence 
ui is indeed analytic in Text- If 00 e Text, we find by a similar argument for the expansion 
of w at 00 



21, Jr aonC)...<-i(C) z^Q 



«o^(c)...a,^_i(c)(^-c); ^ v<(c)-.-«^i(c) 

which again by Theorem l2 .101 starts with the term z~"~^. □ 

3. Convergence results for multipoint Pade approximants 

The aim of this section is to generalize various convergence results for complex Jacobi 
matrices to the setting of linear pencils. 
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3.1. Numerical ranges of linear pencils. It is well known that zeros of formal orthogonal 
polynomials lie in the numerical range of the corresponding tridiagonal operator. More- 
over, the corresponding sequence of Pade approximants converges locally uniformly out- 
side the closure of the numerical range lfT2l Theorem 3.10]. In this section, we generalize 
this machinery to the case of linear pencils and multipoint Pade approximants. 

Let us recall that, for a bounded operator T acting in its numerical range is defined 

by 

e(T) :={(Ty,y),. : 112/11 = 1}CC 

Clearly, 6(T) is a bounded set. By the Hausdorff theorem we have that the spectrum (t(T) 
of T is a subset of the convex set Q{T) (for instance, see [[251 Section 26]). The following 
definition generalizes the concept of numerical ranges to the linear pencil case. 

Definition 3.1 (|l25l). The set 

W{A, B) -.^ {z e C : {{zB ~ A)y, y)p ^ Ofor some y ^ 0} 
is called a numerical range of the linear pencil zB — A. 

The following proposition is immediate from Definition [3T| 
Proposition 3.2. All the zeros ofqn and Pn belong to W{A, B). 

Proof. Let us suppose that ^ is a zero of the polynomial (j„. Thus, according to ( 12.4b . there 
exists an element € C" such that 

(C^[0:n-1] - ^[0:ti-1])% = 0, = 1. 

The latter relation impHes ^ G M^(^[0:n-i] 7 ^[0:n-i]) C W{A,B). Similarly, we have the 
inclusion of the zeros of pj to B). □ 

In general, for the bounded operators A and B, the set B) is neither convex nor 

bounded. However, it turns out that the condition 



O^e(B) (3.1) 



implies a{A, B) C B) Il25l Section 26], as well as the representation 



W{A,B)^{^f^: f^Q) = {^jll: ||/i| = l| , (3.2) 

from which we see the boundedness of B). Condition ( 13.1b implies that B is bound- 

edly invertible, but, in contrary to the spectrum, in general it does not imply any link with 
the numerical range of the operators B^^A or AB^^ . To see this, one may extend the 
following simple 2x2 example to the pencil case 

A 

where e(J5-iA) £ Ql.B-^/'^ AB-^/"^) = <d{A,B) £ e(AB-i) are ellipses with the 
same foci but different eccentricities. 

Generalizing lfT2l Theorem 3.10] for complex Jacobi matrices, we are able to prove 
a result on locally uniform convergence which in some sense generaUzes the Gonchar 
theorem |fT9l . 

Theorem 3.3. Let ( 13.1b be satisfied. Then the sequence of multipoint Pade approximants 
"^[o:rl] Pn+il<ln+i Converges to the m-function locally uniformly in C \ B). 



" 1 1 ■ 


B = 


■ 1 


" 


1 





2 
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Proof. Denote by D C C \ B) a closed set with compact boundary. Setting d 

Jnf /)| > 0, we find for z £ dD and ||/|| = 1 that 

\\{zB^A)f\\>\{BfJ)\ 
implying that 



{AfJ) 



(BfJ) 



> d dist{z,W{A,B)), 



max II (zB- A)- 



< di := — max 



1 



dzedD dist{z,W{A,B))' 



Since M^(A[o:„_i], i?[o:n-i]) C B), the same argument can be used to estimate the 

norm of the resolvent of finite subsections 



max 

z£dD 



(zBu 



5[0:n] - A[0:7l], 

Let i/j be a finite sequence, that is, ip ~ {tpi , . . . , V^fc , 0, 

{zB - A)'tjj = {zB[o,j] - A[o:j])V' = 77 
for sufficiently large j e and 77 is also a finite sequence. Further, one obviously has 



,)^. Then 



(3.3) 



(zB-Ay^r] 



lim {zBio:j] 



^[0:,])"S. 



(3.4) 



Since zB — A is bounded and boundedly invertible, the set of such t^'s is dense in and, 
therefore, due to ( 13.31 ) we have that formula ( 13.4b is also valid for all jy e £^ implying the 
pointwise convergence m[o:j] (z) — s- to(z) for any z G C\W{A, B). Now, the statement of 
the theorem immediately follows from (13.3b and the Vitali theorem Il32l Section 5.21]. □ 



Notice that the concept of a numerical range is valid for operator- valued functions ll25l . 
Thus the presented approach can be also generalized to linear pencils proposed in ||6]. 

3.2. Uniform convergence of subsequences in neighborhoods. We start by improving 
the pointwise convergence result of Corollarv 12.61 generalizing |j3] Corollary 3]. It was 
Ambroladze 12] Corollaries 3 and 4] who first observed that, for real Jacobi matrices, a 
quite dense subsequence of convergents of ( 12.1b converges uniformly in a neighborhood of 
any element of the resolvent set. This result has been generalized in IfTTl Theorem 4.4] to 
the setting of complex Jacobi matrices. We follow here the lines of the proof presented in 
||9] Theorem 4.7] since this allows to deduce in the next subsection a result of convergence 
in capacity in bounded connected components of p{A, B). 

A central observation in what follows is the following result which for complex Jacobi 
matrices may be found in IfTTl Proposition 2.2]. 



Proposition 3.4. The family of rational functions 

1n{z) qn{z) 



u„{z) 



g„+i(z) ali{z)q^^^^{z) a^{z)q^^^^{z) 
is normal with respect to chordal metric on p{A, B). 

Proof. We only have to show that w„ is equicontinuous on the Riemann sphere. By the 
definition of the chordal metric we find for x,y £ p{A, B) 





an{x)qn+i{x)q^iy) - a^{y)qk{x)q^+i{y) 








[<iniy)^aniy)<in+iiy)] 



In order to minorize the denominator, we write shorter as in the proof of Lemma |Z7 
(q^,..., 0, 0, . . . ), = [q^, , 0, 0, ... )^ 



1[0:n] 
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and observe that 



<l[0:n]ix)ixB - A) = [^0^^^,a„(x)g„+i(x),-a„(x)g„(a:),0,. 

n 

implying that 
Similarly, 



{yB - A)<z^;„](y) = 0, 0, a^^(j/)g//+i(y), 0, .. 



T 



implying that 
Finally, 

= 9[0:n]i^)[ixB -A)^ [yB - = - y)qiO:n]ix)Bqfo:n]iy). 

and a combination of these findings yields that x(un(a;), Un{y)) is bounded above by \x—y\ 
times a quantity which can be bounded for x, y lying in compact subsets of p{A, B). □ 

We are now prepared to generalize fTTl Theorem 4.4] to linear pencils. 

Theorem 3.5. For any ^ G p{A^ B) there exists a closed neighborhood V C p{A, B) and 
e„ € {0, 1} such that rn^^Q-n-i+e,^] converges to m uniformly in V. 

Proof. Let f„ = m„ and e„ = if |Mn(OI < !> or elsewhere w„ = 1/un and e„ = 1. Then 

\m{z) ~ TO[0:„-1+£„](2;)| 











(^) 





Using the equicontinuity of the it„ (and thus the w„) established in Proposition 13. 41 there 
exists a neighborhood y of ^ such that li^nC^)! < 2 for all z . Applying the Cauchy- 
Schwarz inequality to ( |2.19t . we obtain for z £V the upper bound 

\m{z) - TO[0:„_1+<:„](2;)| < 

y5y^|r^(z)P + |a^Wr,i+i(z)|y|r«(z)P + \a^{z)r^^,{zW 
and the right-hand side tends to zero with a geometric rate according to Remark l278] □ 

One may construct examples with B = I and selfadjoint A with the spectrum C \ 
p{A, B) consisting of two intervals being symmetric with respect to the origin ^ = 0, and 
m[o:„_i] has a pole at ^ for all odd n. This shows that we may not expect convergence for 
a subsequence denser than that of Theorem l3.5l 



3.3. Convergence in capacity. As explained already before, in general one may not ex- 
pect convergence of TO[o:n] to m locally uniformly in p{A, B) since there might be so- 
called spurious poles in B). One strategy of overcoming the problem of spurious 
poles is to allow for exceptional small sets, as done in ifTTl Theorem 3.1] for complex 
Jacobi matrices where convergence in capacity is established. We may generalize these 
findings for linear pencils, where again we follow the lines of the alternate proof presented 
in iH Theorem 4.7]. 
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Theorem 3.6. Let V be a closed connected subset of p{A, B) with compact boundary, 
then there exist e„ G {0, 1} such that m[o:n-i+e„] converges to m in capacity in V. 

If ( 13.1b is satisfied and V ^ B) then we obtain convergence in capacity of the 

whole subsequence. 

Proof. Let again be u„ = with e„ G {0, 1} to be fixed later, and consider the sets 

V,:^{zeV: \v„{z)\ > 1/e}. 

The arguments in the proof of Theorem 13.51 show that m[o:„_i+e„] converges to m uni- 
formly mV\Ve. It remains thus to show that the capacity of Ve tends to zero for e — > 0. 

We choose e„ in order to insure that the normal family does not have a partial 

limit being equal to the constant oo in the connected component of p{A, B) containing V: 
this can be done for instance by choosing a fixed ^ €V and to take e„ as in Theorem l3.5l 
namely e„ = if |'u„(^)| < 1, and elsewhere e„ = 1. However, under the assumptions of 
the second part of the statement, by taking ^ G F \ VK(y4, B) it follows from the proof of 
Theorem l3.3l that 



sup |u„(OI = sup |e^(^B[o:„] - ^[o:„]) 



< oo. 



and hence here we may take the constant sequence e„ = 0. 

It is now a well-known fact on normal families (see for instance fTT', Lemma 2.4] or the 
proof of 12] Theorem 4.7]) that for normal meromorphic families with partial limits 

different from oo there exist monic polynomials uJn of bounded degree independent of n 
such that 

C := sup max |w„(z)'i;„(z)| < oo. 



This enables us to ensure that 

C 



Ve^{zeV 



> 1/e} C {z G C : |a;„(z)| < eC}. 



|w„(z)| 

Since the capacity increases for increasing sets, and since the capacity of the right-hand 
lemniscate can be explicitly computed to be (eC)^/ dcgLj„^ j-j^g assertion is proved. □ 



4. BlORTHOGONAL RATIONAL FUNCTIONS AND BI-DIAGONAL DECOMPOSITIONS 

In this section we give an operator interpretation of the Darboux transformations of 
rational solutions of the difference equations in question (for the orthogonal polynomials 
case see fl3l ). In other words, we present a scheme for constructing biorthogonal rational 
functions. As a special case, we can construct orthogonal rational functions. Note that 
more information about orthogonal rational functions can be found in llT4l . 



4.1. LU -factorizations. Let us try to factorize the linear pencil zB 

zB- L{z)D{z)U{z) 



A as follows 



(4.1) 



where D{z) ~ diag((io(-z), di{z)^ 
ces of the forms 



..) is a diagonal matrix, and L, U are bidiagonal matri- 



L 



( 1 



V 



U 



( 1 



V '■■ 



Comparing coefficients gives 



do — Po, Pn 



dn- 
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Thus do{z) ~ qi{z)/qo{z) by ( 12.31 ). and by recun-ence using ( 12.2b one deduces that 

, / ^ <ln+l{z) L N _ «n(^)gn(^) R, n _ «n(^)g«(^) 

qn{z) qn+l{z) qn+l{z) 



Hence, the decomposition ( 14.11 ) exists if and only iff qn{z) ^ for all n > 0. In particular, 
from Proposition ^ .21 we obtain existence of such a factorization for z ^ W{A^B). 

The decomposition ( 14. Il l gives us the possibility to define Christoffel type transforma- 
tions. 

Proposition 4.1. Under assumption ( 12.221 ), let xq e Text such that the decomposition 
( 14.71 ) exists for z = xq. Define for n > the functions rational in x 

Q^[xo,x) = ^ , Q^{xo,x) = ^ . 

Xq — X Xq — X 

Then we have the orthogonality relations 



1 



— / QAxq,x)Q^{xq,x){xo - x)m{x)dx = 5j^k/dj{xQ), (4.2) 
where Sj,k is the Kronecker delta. 

Proof. Denote by /j ^ the expression on the left-hand side of ( 14.21 ). We only consider the 
case < j < k, the other case follows by symmetry. By definition of Qj'{xo,x) and 
Qj^{xo, x) and by Lemma lZOl we obtain 

Ijk = R{xo)k,j - Vkixo)R{xo)k+ij 

-vf{xo)R{xo)k,j+i + vf{xo)vk{xo)R{xo)k+i,j+i- 

For j < k, we may apply Theorem l2.3l and obtain after factorization 

Ijk = (q^ixo) - vfixo)q^+,{xo)){r§{xo) - v^{xo)ri+i{xo)) = 

by definition of vf{xo). If j = k, we get slightly different formulas from Theorem l2. Bl and 
obtain after some simplifications 

^ qfix,)irfixo) - (xo).- ,(.o)) ^ ^ = 

where in the second equality we have applied (12.191 ). □ 

Remark 4.2. Clearly, the functions afj . . . a^^Q^^{xo, •) and . . . a^-Q^{xQ, •) are poly- 
nomials of degree < n. □ 



Proposition ^. 1 I tells us that the Christoffel transformation leads to multiplication of the 
biorthogonality measure m{x) by a linear factor (xq — x). This process can be repeated. 
Indeed, after the Christoffel transformation we again obtain a pair of biorthogonal rational 
functions satisfying a generalized eigenvalue equation with a new pair of the Jacobi ma- 
trices A, B 1341 . We can thus apply the Christoffel transformation to these new functions 
factorizing the linear pencil .tiS — ^ in a similar way as in ( 14.11 ). Then the weight function 
m{x){x[) — x) is multipUed by a linear factor xi—x with xi ^ xq. Repeating this process, 
let us introduce the polynomial 7r7v(a;) = [xq — x){xi — x) . . . {x^q^i — x) with Xi ^ xj, 
for i ^ j and construct the functions 

Q„[xo,xi, . . .,xn-i;x) 



■KN{x)Bn,N 



16 



BERNHARD BECKERMANN, MAXIM DEREVYAGIN, AND ALEXEI ZHEDANOV 



where 



<iv(^) = det 



L 

■i+l 



dot 



'+1 



(4.3) 



(4.4) 



and similar expressions for Q^{xo, Xn-i; x), j^{x) and B^j^{x). Note that if two 
or more of the parameters Xi coincide, say xi = xo, then we may apply a simple limiting 
process leading to appearance of derivatives in corresponding determinants. Then it is easy 
to show that these functions satisfy the biorthogonality relation 



. ,XN^i;x)Q^{xo,xi, . . . ,XN~.i\x)iTj{x)m{x)dx 

,Xj-i), (4.5) 
. . , Xj-i). Formulas ( I4.3l )-( l4.5l l are a direct generalization 



= 5j_k/dj{xo,xi,. . . 
with some constants dj (xq, xi, . 

of the Christoffel formula for the orthogonal polynomials, see, e.g., [STl §2.5]. 
4.2. U L-decomposition. For z G p{A, B), let us find a decomposition 

zB - A = U{z)D{z)L{z) (4.6) 
with a diagonal matrix D{z) = diag(rfo(-z), di (z), ■■■■), and bidiagonal matrices 



U = 



/ 1 


-< • 


•• \ 




/ 


1 








•• \ 





1 -uf 










1 
























1 




V : 




■■ J 




V 








■ / 



By comparing coefficients we have 

-a« = -Undn+i, -an = 



-U^dn+1, Pn = dn + U^U^'dn+1 = d,i + 



71 



dn 



+ 1 



It turns out that this decomposition is unique after fixing an arbitrary value for do. Indeed, 
let U-i = do, 2/0 = 1- and consider ?/„ defined by the recurrence relation (I2.2l i. Then it 
follows that 



Vn 



Vn+l 



Vn+l 



where from ( I2.3l l we learn that 

y„(z) = (1 - m{z)dQ{z))qn{z) + do{z)rn{z) = g„(z) - do{z)pn{z) (4.7) 

for all 71 > —1. Thus the decomposition ( 14.61 ) exists if and only if d[){z) g C is chosen 
such that yn{z) 7^ and a^^{z)a^{z) ^ for all n > 0. For the special case d[j{z) = we 
may compare with the LU decomposition of the preceding subsection and get ~ 1 /v^ 
and similarly = Also, for the special case m{z)do{z) = 1 one may show that 

yn{z) = do{z)rn{z) ^ provided that z ^ W{A,B). 

Suppose that the above factorization exists for z ~ xq, and define the Geronimus type 
transformations by the following formulas 



Qnixo,x) ^ q^[x)-u^_^{xo)q^^^{x), Q^{xo,x) 
and Q^{xq,x) = Q^{xa,x) = 1. 
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Proposition 4.3. Under assumption (12.22b . letXQ G T^xt, do{xo) ^ 0, such that the above 
factorization (14.61 1 exists for z = xq. Consider for g G C(r) the linear functional 

&i9) = ^ / 5(0^^ dC + (-^ - m(xo))5(xo), 
then we obtain for j ^ k the biorthogonality relations 

©(Qf(-,xo)Q^(-,a;o)) =0, G(Qf{.,x^)Q^{-,x^))=j^^. (4.8) 

Proof. We only look at the case j > k > 0, the other case follows by symmetry. Let us 
compute the {j, fc)th entry of the product L{xq){xoB — A)^^ (which formally is perhaps 
expected to be equal to the upper triangular matrix D{xq)^^U{xo)~^ but turns out to be a 
full matrix). Using Lemma IZ91 and observing that xq £ Text we get for j > 0, 



{L{xo){xoB ^ A) ^ek,ej) 
= {{xqB - Ay^ek,ej) - uf_-i^{xo){{xQB - A)~^efc, ej_i) 

Note that the same conclusion is true for j = 0. If now j > k, we may rewrite the last 
expression as 

{L{xo){xoB - A)-iefc,ej-) = (rf{xo) - uj^_i(a;o)rf_i(xo))<7^(a;o). 

Noticing that uf^^{xo) = (a;o)/yjii(a;o) with ^ i/„/(a^...a^_i), we get accord- 
ing to (|477] | 

rf{x„) - itjli(.To)rjii(xo) 

■' dQ[xQ) ■> y J do(a;o) y 

Thus for all g £ spanjg^ : fc = 0, j — 1} = span{Q^(-, .tq) : fc = 0, — 1} we 
conclude that &{Qf {-, xo)g) = 0, and, by definition of & and Theorem l2.3l 



{L{xa){xaB - Ay^ej,ej) + ( , / , - m{xo)) Qf{xo, xo)q'f (xq) 

(a;o,xo)(r[(a;o) + (-7-^ — - - m{xa))qf {xo)) 
\ J do(XQ) J J 



Sf(xo,xo)- 



' do{xa)a^{xQ)...a'-_^^{xa) dj{xo)' 
the last claim being evident for j = 0, and for j > according to (I2.6l l and ( 14.71 1 

Qf{xo,xo) ^ pjg.,-_i - Pj-iQj 

do{xo)a^ {xo)...af^j^{xo) yj^ia^ ...a'-_^a^ ...af^^ 

1 _ vAxg) 



yj-iaf^iaf_-^ dj{xo) 
where for simplicity we have dropped in the intermediate expression the argument xq. □ 
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Remark 4.4. Formula (14.81 ) means that the (bi)orthogonality measure m{x) for the trans- 
formed rational functions Qj'{xo, x), Qf{xo, x) consists of a regular part m(x)/{xo — x) 
on r plus a point mass at x = xq, with the mass Mq = 27ri(l/do(a;o) ^ rn{xo)), where 
'^o(2;o) 7^ is a free parameter. A similar situation occurs in the case of ordinary orthogo- 
nal polynomials, where the additional point mass in the Geronimus transformation can be 
freely chosen ITS). □ 

Remark 4.5. Proposition 14. 3 1 for Xq —>■ oo (after multiplication with xq) has been consid- 
ered before in [|18l Theorem 2.2]. □ 



5. An example 



In order to illustrate the above findings and to give a non-trivial example, we study in 
this section the properties of a symmetric linear pencils related to a Markov function of the 
form 



tp{z) 



t ' 



with a probability measure ^ with a support included in some compact real interval [a, h]. 
Here, the entries Aj^k, Bj ^ of the hnear pencil zB — A for symmetric interpolation points 

zi = z^,Z3 = z^, ... e C\[a,b], dist{zj,[a,b]) > 5 > 0, j € N, (5.1) 

are obtained by developing ip into an even part of a Thiele continued fraction. Before 
going into details, we recall from the beginning of SjT]the special case of interpolation at 
infinity zi = Z2 — z^ = ... = oo. Here the expansion of Lp into a J-fraction generates a 
pencil zB — A with a real Jacobi matrix A and with B ~ I the identity, and it is known 
that the spectrum of the linear pencil zB — A (and thus of A) is given by the support of 
the underlying measure /x, and the numerical range equals to its convex hull [a,b]. The 
aim of this section is to show that these properties remain valid for more general sets of 
interpolation points. 

Returning to the task of developing ip into the continued fraction in question, the fol- 
lowing result has been shown in IfTSl Lemma 3.1 and Remark 3.3], by making the link with 
Nevalinna functions. The proof given in ITSl uses the assumption | Im Zj\ > 6 > and it 
can be immediately generalized to our setting. 

Proposition 5.1. Suppose that ( 15.11 ) holds, and that ^ has an infinite number of points 
of increase such that p is not a rational function. Then there exist probability measures 
/iQ = M, Ml , Ai2 , • • • such that, for all j > 0, 

tpj (z) = 



^Bjj — Ajj B'^_ 



j+i.j" 



(z - Z2j + l){z - Z2j+2)^] + l{z) 



with the Markov functions 

z-t'- 



LP J [z) 



and the real numbers 



Bj,] - 



f 

J a 



b dfij (t) 



Z2j + l-t 



> 1 



rb tdfj,j(t) 
Ja \z2j + i-t\'^ 



J a 



b dfij{t) 



Z2j + l-t 



B 



j+i 



, = VBJJ^ > 0. 



Hence, our Markov function ip for the symmetric interpolation points ( 15.1b induces a 
linear tridiagonal pencil zB — A if we set according to ( 12.9b . 



zB 



A 



af{z) = zBj+i^j - Aj+ij = - Z2j+i) 



-af{z) 



^BjJ+l — ^jj-+l — Bj + ij{z — Z2j+l) — Bjj^i{z — Z2j+2)- 
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We collect some elementary properties of this pencil in the following two propositions. 

Proposition 5.2. Suppose that ( 15.1b holds. Then the above tridiagonal matrices A, B are 
Hermitian and bounded. 

Proof. It follows from ( 15.11 ) and the explicit formulas given in Proposition 15 . 1 1 that A and 
B are hermitian, and B is real. In order to show that B is bounded, it is sufficient to show 
that its entries are uniformly bounded, where in our case it is sufficient to consider the 
diagonal ones. Let us first establish the minorization 



t 



> 



dist(z, [a, b\ 



z eC \ [a,b]. 



(5.2) 



max{|z — aP, \z — 6p} ' 

For a proof of ( 15.21 ) we suppose that Rez > {a + b)/2, the other case is similar. Since 

1 1-^ Im(l/(z — <)) does not change sign on [a, b], we get 



Im 



Im ■ 



llmzl 



Hence our claim i5.2\ follows provided that | Imz| = dist(2, [a, b]). Otherwise, we have 
that Re(2 - t) > Rc{z - b) > for all t € [a, 6], and hence 

and the claim follows by observing that |z — 6| = dist(z, [a, b]). 

Combining (15.2b with the definition of Bjj given in Proposition ^. ll we conclude that 

^ ^ max(|z2j+i - a\^, \z2j+i - ^l"^) 
^'^ - dist(z2j+i,[a,6])4 

the right-hand side being bounded according to assumption (15.1b . Thus B is bounded. 

Similarly, one shows that the diagonal entries Ajj of A are uniformly bounded. In 
order to discuss the off-diagonal entries of A, we choose a fixed point z G C\[a,b] having 
a positive distance from the set of the interpolation points Zj, and get with the help of 
Proposition 15. II 



1^, 



l^w+iP 



\Z23A 



,2^2 



< 



1 



+ 11 



zB, 



A, 



1 



iPj + l{z) |z-Z2j + l|2 

the right-hand side being bounded uniformly for j > according to (15.2b . Hence, A is 
also bounded. □ 

Proposition 5.3. Suppose that (15.1b holds. Then for all y — (yo, ...)^ e there holds 
{By,y)>\yk\^ ifyo = ... = yk-i = 0. (5.3) 
Furthermore, for the numerical range of Definition 13. 1 1 there holds W{A, B) C [a, b]. 
Proof. In order to show ( 15.3b . let y ~ {yg, 2/1,2/2, --O^ G i'^. We write as before j/[o:„] = 



(2/o,yi,...,2/„,0,0, ...) ' e ^jjjj jjQtjj.g that (By[o 



since B is bounded by Proposition l5.2l Then for yo 
relation Bj j = 1 + j we get 



{By,y) for 



yk-i = and n> k using the 



{By[0:n],y[0:n] 



BjM? + 2 ^j+i^j- i^c(%%+i) 



j=k 



\yk\ 



ri-1 



\Bj+i,jyj + 2/j+il' + Bl+ijyn\^ > \yk\ 

J=k 
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implying that (15.3b holds. Since also (A2/[o:ri]: y[0:ri]) ^ {^yiV) for n — > oo, we get 
using (15.3b and (13.2b that B) is included in the closure of the union of the numerical 

ranges VF(A[o:„], S[o:ri]) of all finite sections. Further, observe that 

W^(^[0:„],S[0:„]) = 0(i?[O:t]^[O:n]S[oi]), 
_ 1 _ 1 

where -SjQ.^jyl[o:„]i3jQ.^j is self-adjoint. Thus, M^(A[o:„] , i3[o:ri] ) is a convex hull of eigen- 
values of the matrix _BjQ.^jyl[o:„]i?jQ.Jjj or, equivalently, the zeros of the polynomial 

dct(z-Bjp^^jA[0:„]SjQ^^j) = detSjQ^^^j dct(zB[0:„]-A[0:„]) = dctBj„^^jq„+i(z). 

To complete the proof it remains to note that all n + 1 roots of qn+i, that is, the poles of 
a rational interpolant of a Markov function are lying in the open interval (a, b), see ll30l 
Lemma 6.1.2]. 

□ 

The positive definiteness of finite sections of B also for not necessarily bounded [a, b] 
has been shown already in ifTSl Proposition 4.2], where the authors also establish (15.3b . 

Notice that property (15.3b in general does not imply that condition (13.1b is true. How- 
ever, only the latter condition allows us to conclude that the spectrum of the pencil zB — A 
is included in [a, b]. There is a special case where we may say more. 

Theorem 5.4. Beside (15.1b . suppose in addition that Z2j+i = ~> 00 fl^^ j ~^ 00. 

Then the operator B is a compact perturbation of the identity, and condition ( 13.1b holds. 

In particular, the spectrum of zB — A is given by the support of the measure /i, and, 
outside the spectrum, ip coincides with the m-function of the linear pencil zB — A. 

Proof. We have shown in the proof of Proposition 15.21 that lAj+ijp = \z2j+i\'^(Bjj — 
1) = \z2i+i P^j+i j is bounded for j 00, and hence 

lim _Bj-|_i j = lim Bj j^i = 0, lim Bj j = 1, 

showing that S is a compact perturbation of the identity, and B has its numerical range 
included in [0, +00) by ( 15.3b . Hence, if ( 13.11 ) does not hold, then would be an eigenvalue 
of B, with corresponding eigenvector y E i'^ ,y ^ 0. Inserting this y into ( 15.31 ) with k such 
that yk gives a contradiction. 

It follows from the text after ( 13.11 ) together with Proposition 15.31 that the spectrum 
a{A,B) of the linear pencil zB — A is included in [a,b]. Also, by construction and 
Corollarv l2.12l Pn/qn interpolates both ip and m in Z2n, implying that these functions 
are equal for z = Z2„ and for all n, and analytic in C \ [a, b] including 00. Since these 
points accumulate at 00, we conclude that m = (p outside [a,b]. Finally, the inclusion 
supp(/i) C (t{A, B) follows from the fact that ip is not analytic in any domain containing 
points of the support of p. 

Given z G C \ [a, 6] , by choosing a contour T surrounding [a, b] but not the interpolation 
points Zj nor z we get from Lemma l279l the formula 

{{zB A)-W,e,) = ^ |^,«(C),,^(C)^dC ^ l\fit)q^,it)^, 

where for the second identity we have used the Fubini theorem and the fact that ip ~ m 
on r. We denote by R{z) the infinite matrix with entries R{z)j^k given by the above 
right-hand integral, which is clearly well defined for any z outside the support of p. From 
Lemma 122] we know that R{z) is a formal left and right inverse of {zB — A), and the 
desired conclusion z ^ cr{A, B) follows as in the proof of Theorem l2.3l bv showing that 
R{z) is bounded. 
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For this last step, we consider the U L decomposition of B discussed in Remark l43] and 
in ifTSl Theorem 2.2]: let U be an upper bidiagonal matrix with ones on the diagonal, and 
the quantities Bj^i j on the main upper diagonal, then U represents a bounded operator 
on P according to Proposition 15. 21 Moreover, we have that B — UU* , and, with B, also 
U has a bounded inverse. Hence it will be sufficient to show that 



\{U*R{z)Uv,y)\ < 



{y,y) 



(5.4) 



dist(2:, supp(//)) 

for all y = {yo, yi, y„, 0, 0, ...)^ G and for all n. Comparing with Proposition 14.3 
we find that 

{U*Riz)Ue,,e,) = f Qf[^,t)Q]:{^,t)'^ 

J a Z l 

where 

Q^,{go,x) = g^(a;) - B„,„_ig,t_i(x), Q^^{oo,x) ^ q^{x) - B„,„_i9^_i(a;), 
and Qq(oo,x) = Qq{oo,x) = 1, and finally 
1 

In addition, since (/„ has real coefficients, it also follows from (15.1b that, for t E 



implying that 



\{U*Riz)Uy,y)\ 



J2yjQf{^,t) 



3=0 



dfi{t) 



< 



1 



dist(z, supp(/i)) 
as claimed in ( I5.4l i. 



{y.y) 



dist(z, supp(/^)) ' 



□ 



Remark 5.5. The assumption Z2j+i = z2j+2 - 
be relaxed. For instance, if 

^{\z2J+l -a|,|z2j+i -^1) 



lim sup ■ 



dist(z2i+i, [a, 6]) 



oo as j — > oo is very restrictive and can 



< </2, 



then it follows from the proof of Proposition l5.2l that sup Bj+ij < 1. As a consequence. 



the operator U from the proof of Theoi'em l5.4l and thus i? is a boundedly invertible operator. 
This implies that (13. Il l and hence the second part of the statement of Theoi'em l5.4l is still 
true. □ 



In the setting of Theorem 15.41 we may therefore apply our findings of Theorem 13.31 
Theorem 13.51 or Theorem 13.61 in order to study the convergence of the multipoint Pade 
approximants towards the Markov function Lp, compare with ifTSl Theorem 6.2]. 

Finally, returning to the discussion of Remark lTTj concerning the degrees of freedom of 
representing multipoint Pade approximants via linear pencils, it is not difficult to see that 
the two linear pencils zB — A and AD{zB — A)D^^A for diagonal D, A with non-zero 
diagonal entries generate the same continued fraction (12.1b . Notice that the matrix D does 
not affect the diagonal entries and can be therefore be considered as to be a balancing factor 
for the offdiagonal entries, whereas A allows to scale the entries. In terms of the continued 
fraction (12.1b . a scaling corresponds to considering an equivalence transformation of ( 12.1b . 
and different normalizations can be found in the literature concerning the special cases 
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of J-fractions, T-fractions or Thiele continued fractions. A balancing, however, leaves 
invariant the continued fraction (12. Il l and just addresses the question how to factorize the 
products oi'-a^. 

It is always possible to choose a scaling such that the resulting matrices A, B become 
bounded. However, such a scaling might produce a matrix B having no longer a bounded 
inverse, or satisfying no longer the condition ( 13.11 ). We also know from fT?, Theorem 2.3] 
that, for fixed z, the balancing which is best for obtaining z e B) is the one which 
makes zB — A to be complex symmetric (i.e., a complex Jacobi matrix). In the special 
case of Theorem 15 .41 we have chosen a balancing factor to make B real symmetric, and a 
scaling such that A, B are bounded and B has a bounded inverse. 

A study of best scaUng or balancing for general linear pencils is beyond the scope of 
this paper For future research it might be interesting to consider a (formal) factorization 
zqB — a = Mi{zo)M2{zo) for some fixed zq (as done in § 4) and to discuss the conver- 
gence of multi-point approximants in terms of spectral properties of z !-> Mi{zo)^^{zB — 
A)M2{zo)~^ , since this latter quantity does not depend on scaling or balancing (but de- 
pends on how to choose the factors Mj (zq)). 
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